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Abstract 

We study Hartree-Fock, Gutzwiller, Baeriswyl, and combined Gutzwiller- 

Baeriswyl wave functions for the exactly solvable one-dimensional 1/r- 

Hubbard model. We find that none of these variational wave functions is 

able to correctly reproduce the physics of the metal-to-insulator transition 

which occurs in the model for half-filled bands when the interaction strength 

equals the bandwidth. The many-particle problem to calculate the variational 

ground state energy for the Baeriswyl and combined Gutzwiller-Baeriswyl 

wave function is exactly solved for the 1/r-Hubbard model. The latter wave 

function becomes exact both for small and large interaction strength, but it 

incorrectly predicts the metal-to-insulator transition to happen at infinitely 

strong interactions. We conclude that neither Hartree-Fock nor Jastrow-type 

wave functions yield reliable predictions on zero temperature phase transitions 

in low-dimensional, i.e., charge-spin separated systems. 
PACS1993: 71.27.-ha, 71.30.-hh 
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I. INTRODUCTION 



Variational wave functions are often used to study ground state properties of quantum 
mechanical many-particle systems. Examples of approximate ground state wave functions 
for correlated Fermions are Hartree-Fock [|l| , Gutzwiller p| , local Ansatz , and, in general, 
Jastrow-Feenberg wave functions [Q. One of the advantages of variational wave functions 
is the fact that they give an exact upper bound for the true ground state energy. Thus 
there is, at least in principle, a criterion to assess variational wave functions: the lower the 
variational ground state energy, i.e, the tighter the variational bound, the better the wave 
function ("energy criterion"). 

In practise, however, two major problems arise: (i) the energy might not be a good 
criterion for the "quality" of a wave function, i.e., a wave function might give a very good 
approximation for the ground state energy but it may still miss the ground state physics 
of the system under consideration. This issue can, however, only be addressed, if the exact 
solution of the corresponding Hamiltonian is known; (ii) the evaluation of correlated wave 
functions poses yet another quantum many-particle problem such that even seemingly simple 
wave functions (and their straightforward improvements) are not analytically tractable. 

To deal with these two problems we focus on the exactly solvable one-dimensional 1/r- 
Hubbard model [^]. It describes spin-1/2 Fermions on a chain of L sites, hopping with long- 
range amplitude, ti^rn = it{—lY~'^ [d{l — m)]"^ = t^^. Here, d{l—m) = {L/n) sin[7r(/— m)/L] 
is the chord distance between sites / and m on the chain closed into a ring (the lattice 
spacing a is set to unity). The hopping becomes ti^m —>■ it/{l — m) in the thermodynamical 
limit L ^ oo for fixed distance (/ — m) ("1/r- hopping"). The electrons interact through a 
local Hubbard interaction, U, and the total 1 / r-Hubbard Hamiltonian reads 

H = f + UD= krn&t^rn,. + uj2ni^rni,l ■ W 

l^m=l,a- 1=1 

For even L we choose antiperiodic boundary conditions, so that the resulting dispersion 
relation is linear in wave vector, namely, e{k) = tk with k = A(m -|- 1/2) (A = 27r/L, 
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m = —L/2, . . . , L/2 — 1). For U = the Fermi sea is the ground state with all fc-states from 
k = — TT to kp = 7T{n~ 1) filled where n = {N^ + Ni) /L is the total particle density. For half- 
filling [n = 1) and U oo, the Hamiltonian ([^) reduces to the spin-1/2 (l/r)^-Heisenberg 
or Haldane-Shastry model |^. In the rest of the paper we will restrict ourselves to n < 1. 
The more than half-filled case can be obtained by particle- hole symmetry 0. 

The exact solution of the 1/r-Hubbard Hamiltonian was conjectured in Ref. 0, and 
its physical properties were discussed in detail elsewhere [^. The model displays a Mott- 
Hubbard metal-to-insulator transition at half-filling when the interaction strength, U, equals 
the bandwidth, W = 27rt. This is in contrast to the one-dimensional Hubbard model which 
was exactly solved by Lieb and Wu |^. There, the tight-binding cosine dispersion has a 
perfect nesting wave vector q = n such that eth{k) = —2tcosk = — etb(A; + g). Consequently, 
the Hubbard model describes an insulating state for all f/ > at half-filling. This metal- 
to-insulator transition at f/ = O"*" for n = 1 is trivially reproduced in an antiferromagnetic 
Hartree-Fock theory [|l|. We therefore use the 1/r-Hubbard model to assess the quality of 
the Hartree-Fock approach for models without the perfect nesting property. 

Another advantage of the 1/r-Hubbard model comes from the fact that we are able 
to calculate variational ground state energies for correlated wave functions without further 
approximations. Thus far, this was only achieved for the Gutzwiller wave function in one 
dimension |]lO|-[T^ , or, for general Gutzwiller-correlated wave functions, in the limit of large 
dimensions [0-|l5| or large orbital degeneracies [Q. For the 1/r-Hubbard model we are 
able to derive the ground state energies for the Baeriswyl [T^ and a combined Gutzwiller- 
Baeriswyl wave function. 

The plan of the paper is as follows: in section || we diagonalize the antiferromagnetic 
Hartree-Fock Hamiltonian to obtain the Hartree-Fock ground state wave function and en- 



ergy. In section |ITT| we introduce and evaluate the Gutzwiller-, Baeriswyl-, and Gutzwiller- 
Baeriswyl wave functions. We draw our conclusions on the quality of these variational 
approaches in section |rV|. 
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II. HARTREE FOCK APPROXIMATION 

As usual m we factorize the interaction part of the Hamiltonian (|l|) as 

and allow for a sublattice magnetization m = \{ni^^) — {ni^i)\, alternating in sign on the 
odd/even lattice points of the A-B lattice. The particle density n = (ni^i) + (ni^i) remains 
uniform. Thus we may write 

(V) = ^ + (-l)Vy . (2b) 

We introduce the new quasiparticle operators ak,a, i>k,a within the magnetic Brillouin zone 
-TT < A; < by 

ak,a = cos 4>k,aCk,a + SIB. (f)k,a Ck+7T,<T 

hk,a = - sin Ck,a + COS (j)k,a Ck+n,^ (3a) 

and the inverse transformation reads 

Ck,cT = COS (f)k,aak,a " slu 0fc,<7 &fc,<7 
Ck+7T,a = sin (f)k,a ak,a + COS (f)k,a bk,a ■ (3b) 

The band width of our model is = 27rt, and the transformation angle (j)k,a fulfills 

2Um 

tan2(f)k,a = (^—^ (4a) 

n 1 / W \ 

sin" (f)k,a = o 1 - ^ / — ■ (4b) 

2 \ ^W^ + {2Umf ) 

The operators are constructed such that al^ (bta) creates particles mainly on the A-lattice 
(i?-lattice) sites. The Hartree-Fock Hamiltonian then reads 

-7r<fc<0 

with 



= tk + j(w + U + Um^) ± j^W^ + AU^m? 



(5b) 



The Hartree-Fock charge gap is given by A/z^^ = E^^_^ — E^^q = y/W^ + AU'^m? — W /2. 
The corresponding ground state wave function reads 



(6) 



Finally, we obtain the particle densities for all temperatures T = as 

, 1 



(n 



k,y 



exp {j3{Ej, - /^hf)) + 1 
1 



fMK) = (nil) 



exp [(3{E+ - /ihf)) + 1 
where the Hartree-Fock chemical potential /xhf is obtained from the particle number 
1 2 /"O r _ 

^ = J E = ^ J_ fpDiEk ) + fpDiE^k 

— 7r<fc<7r,(T ^ 

The self-consistency condition reads 



(7a) 
(7b) 



^a) 



m = \ {n^) - {ni,i)\ 
AUm 



VW^+ALP 



27r 



dk 



fFT){Ek ) - hD{El 



(8b) 



This equation gives the magnetization as a function of temperature and interaction strength. 

We note that the Hartree-Fock Hamiltonian (pa)) will give a finite magnetization at 
finite temperatures, if m > at T = 0. However, a phase transition from a paramagnetic 
to an antiferromagnetic phase at some finite Neel temperature does not happen in a one- 
dimensional system [l^. This is one of the well-known shortcomings of the Hartree-Fock 
Mean-Field approach to low- dimensional systems. Furthermore, it gives a T = transition 
to a magnetic phase at a critical value U = U^^{n) for all band-fillings because a finite 
magnetization is the only way to obtain a lower energy than the paramagnetic Hartree-Fock 
phase for large U/W. We thus restrict ourselves to the half-filled case {n = 1) where the 
Hartree-Fock approach gives the best results. For n < 1 the antiferromagnetic Hartree-Fock 
theory certainly fails to give the correct physical picture. 



In the case of half-filling the Hartree-Fock chemical potential follows from eq. (^) to be 
temperature independent. It is given by /ihf = U {1 + m^) /4 which lies in the middle of the 



gap between the upper and lower Hartree-Fock band. The self-consistency equation (|8ED for 
the magnetization does not considerably simplify but stays an implicit equation which has to 
be solved numerically for T > 0. Here, we are only interested in ground state properties. For 
(3 —>■ oo the lower Hartree-Fock band is completely filled {k^^ = 0), and it follows from 
that either m = (paramagnetic Hartree-Fock solution) or 



2U = VW^ + AU'^m? . (9) 
Hence, Hartree-Fock predicts a transition from a paramagnetic metal to an antiferromagnetic 



insulator at Uf^{n = 1) = W/2, above which we have m = yl — (f/^^/f/)^ with the 

typical mean-field exponent of one-half near the transition. The charge gap is A/i[^^ = 

(VW^TWh^ -W) /2 = U - Uf^. The Hartree-Fock gap grows linearly function 

of f/ — U^^ in agreement with the exact solution |^,^ where, however, the true transition 

is at twice the Hartree-Fock value, Uc = 2U^^. It is amusing to note that the density of 

states for spinon excitations develops a van-Hove singularity (formation of the upper and 

lower Hubbard band) precisely aX U = W/2 p|. It is again seen that antiferromagnetic 

Hartree-Fock is sensitive to spin-correlations only. 

The Hartree-Fock variational ground state energy density at half-filling is given by Cg 

(1/^) E-n<k<oEk which gives 

U-W 



.HF 



4 

W^ 



for U < W/2 (m = 0) 



(10) 



for U > W/2 {m = Jl- {W/2Uf) 



16U 

which is compared in fig. |T] to the exact ground state energy 

(U-W f/2 



eo 



4 12W 
W^ 



ioT U <W 
hi U>W 



12U 

It is seen that Hartree-Fock, as an effective single-particle theory, fails to give any contri- 
bution to the correlation energy of order U^/W. There are no non-analytic contributions 
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of the order exp{—W/U) for small U/W as expected for a Hamiltonian without the perfect 
nesting property. 

At large U/W it gives the correct analytical behavior with a smaller prefactor. However, 
it also predicts a finite sublattice magnetization which is, of course, zero in the exact solution. 
Furthermore, the Hartree-Fock gap is temperature dependent and one gets a transition 
from an antiferromagnetic insulator to a paramagnetic metal as function of temperature, if 
U > U^^. The gap in the exact solution is temperature independent, and no transitions occur 
at finite T. Finally, the zero temperature transition itself in Hartree-Fock is driven by spin- 
correlations. This can be seen from the fact that one has a finite sublattice magnetization 
after the transition, and that the transition happens at all fillings. This is in contrast to the 
exact solution where the transition happens only at half-filling, and no magnetic ordering 
occurs. Hence, the transition in the 1/r-Hubbard model is driven by charge-correlations, 
just opposite to the physics of the Hartree-Fock approximation. 



III. CORRELATED WAVE FUNCTIONS 

We now focus on the correlated wave functions whose energy we want to evaluate. Our 
starting point is the Gutzwiller wave function [Q] 

L 

li^cig)) = S''^ I Fermi-sea) = H [l ~ (1 ~ |Fermi-sea) (12) 

1=1 

which starts from the Fermi-sea, the exact ground state for U = {g = 1). For U > 
double occupancies are less favorite because of the on-site repulsion in the Hamiltonian, 
eq. (p. To include this effect, double occupancies are suppressed by the Gutzwiller corre- 
lator ^ — — g)Di which globally reduces configurations with double occupancies in 
the Fermi-sea, if g < 1. Note that this correlator is regular for all < g < 1, i.e., lipcio)) 
describes a (correlated) metal as long as we do not project out all double occupancies. For 
g = we alter the nature of the Fermi-sea because we not only suppress configurations 
with double occupancies but truly eliminate them with the help of the Gutzwiller projec- 
tor Pz5=o = nf=Ji-A". 



A metal-to-insulator transition within this wave function can only be expected at half- 
filhng, if the variational procedure gives g = such that every site is singly occupied 
(for less than half-filling the remaining holes are still mobile and keep the state metallic). 
Brinkman and Rice observed that such a transition ( "Brinkman-Rice transition" p!9|p0| ) at a 
finite U = U^^ = 8\ Cq^I was contained in the Gutzwiller Approximation [Q] to the Gutzwiller 
wave function, where Cq^ is the kinetic energy of the Fermi-sea in the non-interacting sys- 
tem. It was later understood that this approximation becomes exact in the limit of infinite 
dimensions [|T3,|T3 . If we plainly apply the Gutzwiller Approximation to our one-dimensional 
system we obtain Cq^ = —W/A, Uf^ = 2W, and the ground state energy density at half- 
filling in the Gutzwiller Approximation is given by e^"^ = eQ^(l — U/U^^Y for U < Uf^, 
e'^^ = for f/ > f/^^. This energy density is plotted in figure ^ in comparison to the 
Hartree-Fock and the exact result. Although the Gutzwiller Approximation is an uncon- 
trolled approximation in any finite dimension, we obtain the interesting result that it predicts 
a metal-to-insulator transition at half-filling for a finite value U = Uf^ = 2W. This transi- 
tion is driven by charge-correlations only which is obvious from the fact that the ground state 
energy density in this approximation is zero beyond Uf^. It is thus seen that this approxi- 
mation very accurately describes the physical mechanism of the metal-to-insulator transition 
in our one-dimensional model, despite the fact that it is not variationally controlled, and 
the variational estimate for the ground state energy density is poor for large U /W. In any 
finite dimension, the Brinkman-Rice transition will be shifted to = oo where double 
occupancies are strictly forbidden by the Hubbard interaction. This has indeed been proven 
for the Gutzwiller wave function ||2l[ . 



It was only recently that one was able to go beyond the Gutzwiller Approximation to ex- 
actly evaluate the Gutzwiller wave function in one dimension |[10|JTT[] . It was shown that the 



Gutzwiller wave function does not give a good ground state energy for the one-dimensional 
Hubbard model for large U/W due to poor correlations between double occupancies and 
holes [|T2|. On the other hand, the correlations between spins are excellent such that the 
Gutzwiller projected Fermi-sea is a very good trial state for the one-dimensional antiferro- 
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magnetic Heisenberg model p2l|TI],|T2[ . We will see below that these considerations remain 



valid for the 1/r-Hubbard model, eq. (|lD. 



Baeriswyl [T^] constructed a variational wave function for the one-dimensional Hubbard 
model to overcome the difficulties of the Gutzwiller wave function for large U/W. He started 
from the Gutzwiller-projected Fermi-sea to incorporate its excellent spin correlations, and 
proposed 

= fe^/'^Po^o I Fermi-sea) (13) 

as variational ground state {0 < b < 1). For large U/W the variational parameter b is close 
to unity, and double occupancies and holes are linked by the application of the kinetic energy 
operator to the Gutzwiller-projected Fermi-sea. Hence, the missing charge- correlations in 



the Gutzwiller wave function are now properly taken into account |17 . 

At half-filling, the Baeriswyl wave function describes an insulator for all < 6 < 1 
because the Gutzwiller-projected half-filled Fermi-sea is insulating, and the Baeriswyl- 
correlator 6^/^ is regular as long as 6 > 0. A metal-to- insulator transition will only happen, 
if 6 = is the result of the minimization procedure. Since 6 = corresponds to the free 
Fermi-sea we can expect that the transition will happen at f/ = 0^ for the Baeriswyl wave 
function. Unfortunately, not much is known about the Baeriswyl wave function because 
even the calculation of the variational ground state energy is too difficult for general Hamil- 
tonians. Below, we show how this can be accomplished for the 1/r-Hubbard model. Our 
analysis will confirm our general considerations on the metal-to-insulator transition in this 
wave function. 

Obviously, it is necessary to go beyond Gutzwiller and Baeriswyl wave functions because 
neither of them is able to describe the metal-to- insulator transition of the 1/r-Hubbard 
model at U = W. Their natural generalization is the "Gutzwiller-Baeriswyl" wave function 
defined as 

|7/>gb(6,^)) = &^/^^^|Fermi-sea) . (14) 
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For 6=1 this wave function reduces to the Gutzwiller wave function (0), and to the 
Baeriswyl wave function (|13D for g = 0. We now have a two-parameter wave function 
which allows for the competition of the "kinetic" Baeriswyl correlator and the "potential" 
Gutzwiller correlator. We might hope that this competition will result in a variational 
prediction of a metal-to-insulator at some finite critical value U^^. We will see below that 
this is not the case. Instead, we again find U^^ = oo for the 1 / r- Hubbard model. 

We denote (O) = (^/'gb|0|-?/'gb)/(V'gb|^gb)- Details of the calculation are presented in 
appendix The final result for the expectation value of the kinetic energy is 

and the mean double occupancy is given by 

_ 1 /-sn-i (a(x)\^l + X + B(x)\^l — x] 



where 

a{x) = a{x] b, g) 



9 — 1 / 
1 + ^(1 + ^) 



(16a) 



Vh{g-l 



P{x)^P{x-h,g) = ^^^ '-Vl^ . (16b) 

Before we proceed with the general case, we first discuss the special cases 6=1 (Gutzwiller 
wave function) and g = Q (Baeriswyl wave function). 



A. Gutzwiller wave function 



We set 6 = 1 in eqs. (|T5]), (|^), and obtain a{x; g) = \i — x + g{l + x)] /2, and l3{x] g) 



{g — 1)a/1 — After some straightforward manipulations one arrives at 



4 2 Wo 'l + {g^-l)y 
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A simple integration gives 



^(^) = on ^'2^2 [-(1 - - In (1 - (1 - g^)n)] (18a 



2(1-^7^ 



;i8b) 



4 V^? + iy 

As shown in appendix these results completely agree with those obtained from a direct 
application of the methods developed in Refs. [|l3,|ll|], and provides an independent check 
for the conjectured solution of the 1/r-Hubbard model in Ref. 0. 

To obtain the variational ground state energy density one has to minimize e'^{g) = 
{T{g))/L + Ud[g) with respect to g. This has to be done numerically for general U/W. 
The resulting curve at half-filling is shown in figure It is seen that the Gutzwiller wave 
function is very good for values of U/W < 0.5. Indeed, as shown in appendix ^, the 
Gutzwiller wave function gives the exact ground state energy to order U"^ /W for all fillings. 
It is further exact to order for all values of U /W , see appendix ^. The logarithmic 
dependence of d{g) on g for small g (i.e., large U /W) causes the variational energy to behave 
like Cq ~ W"^/ [U\n{U/W)], see Ref. |l^, in contrast to the exact result, see eq. (pT|). As 
discussed in detail in Ref. [jl2|, this wrong behavior of the variational ground state energy 
comes from the poor correlations between double occupancies and holes in the Gutzwiller 
wave function for large U/W. 

It is clearly seen that (7 > for all f/ < 00 such that the variational prediction for 
the metal-to-insulator transition is U^ = 00. It should be clear, however, that at least 
the physics of this transition is correct: it is due to charge correlations, magnetic ordering 
neither occurs nor is it necessary. Furthermore, the transition only happens at n = 1 as it 
should, and not for arbitrary electron density as in the case of Hartree-Fock theory. 



B. Baeriswyl wave function 



We know that the Gutzwiller-projected Fermi-sea is the exact ground state of the 1/r-tJ 
model with pair- hopping terms i.e., of the Hamiltonian which one gets from eq. ([1|) by 
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first order perturbation theory in W/U. Consequently, the Baeriswyl wave function is the 
exact ground state of the 1/r-Hubbard model (|I]) to order W/U [T^JT^ . 

For g = the equations (p!6D can be simplified to a{x; 6) = (1 - x)/{2Vb), and b) = 
-y/bVT^^/2. After some straightforward manipulations eqs. (^) finally give 

{fm/L = - ^"^^""^ + W—^^ - b')n - In (l - (1 - b^)n)] (19a) 



d{b) 



1 + 6 



2(1-62)2 
{f)/L n{l-n) 



W 



(19b) 



Comparing eqs. (|18a|) and ( |19a| ), and eqs. (|18b| ) and ( [L9b|) , we note that the mean kinetic en- 
ergy and the mean double occupancy just change their roles when we go from the Gutzwiller 
wave function to the Baeriswyl wave function. However, this only happens in the highly 
symmetric 1/r-Hubbard model. 

Form this equivalence it immediately follows that there is a metal-to-insulator transition 
for n = 1 in the Baeriswyl wave function at = 0~^, i.e., this wave function is insulating at 
half-filhng for all U > 0. While for large U/W we obtain the exact result = —W'^/{12U) + 
0{W^ /U"^), the wave function gives only a poor estimate for the ground state energy density 
at small U /W. One obtains the Hartree contribution to first order, but there is, at half- 
filling, no contribution to second order in {U/Wy. Instead, the correction is proportional 
to (UyW) \n{U/W). 

At least for the 1 / r-Hubbard model we can thus conclude that it is as insufficient as the 
Gutzwiller wave function in the opposite regime of its obvious applicability. This unsatisfac- 
tory behavior might in part be due to the special model. Unfortunately, this wave function 
cannot be evaluated for other Hamiltonians away from 6^1. 



C. Gutzwiller-Baeriswyl wave function 

For finite b and g, the Gutzwiller-Baeriswyl wave function represents a metallic state. 
We first analyze whether we can have a metal-insulator transition at a finite value of U /W. 
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The integrals in eqs. (0) can be done analytically. To properly take care of possible 
singularities we write 



{mg))/L 



Wn{2 - n) 



W 



(20a) 



4 2 

2d{b, g) = [l + b{g - l)f (/(I) - + {g- 1)2(1 _ 6)2 (/(3) _ J{2)) + g2j(2) ^20b) 



where 



j(m) 



dy- 



(21) 



A(6,(7)t/2 + 5(6,(?)|/+l 

with A{b,g) = {g - 1)^(1 - 6^) > 0, 5(6,^) = 2{g - 1) + ^^(^ - 1)^ < 0. The denominator 
of the integrals I^"^^ has zeros at Zi^2 = {\B\ ± y/B"^ — AA) / {2A). One can easily show that 
Zi > Z2 > 1- This implies that there will be a singularity (metal-insulator transition) only 
at half-filling (the upper limit of integration is < 1), if Z2 = 1- This can only be fulfilled, 
if A{b, g) + B{b, (?) + 1 = 0, or (7 = 0, in agreement with our general considerations. 

We have to find out whether the variational procedure gives g = Q aX some finite . 
To this end we have to carefully analyze the limit (7 — > 0. We can assume that b will not be 
unity because this would correspond to the case Uf^ = 00. In the limit g —>■ we first note 
that Z2 ^ 1 + 0{g^), ^ 1/(1 - 62). Then, the combinations /(^^ - /(2) and J^^) - /(2) 
stay finite, and actually are of the order g^ In (7. Also, I^'^\b, g) \ng for small g. If we now 
expand eQ^{b,g) around g = 0, one finds 



e™(6, g) = e^ib) + 61(6)^7 + 63(6)^7' In ^7 + 



(22) 



The first term in this expansion, ei(6) = de^^ib, g)/dg ^ exists because all singularities at 
least of the order g^lng. After some lengthy calculations we find 



eiib) 



W 



;i-62)s 



l-b^ + 2bHnb +U 



1 - 52 - 261n6' 



(23) 



Close to the transition we can minimize eq. (p^) with respect to g which gives 

eiib) 



g^T^g 



262 (6) 



(24) 
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Thus we know g{b) as a function of b. We insert this into our expansion (^) and obtain 



e^''ib) = e^{b) + ^-^e^{b) (25) 

For given U /W we have to compare the minima of e[f^(6) at b^^ and of e^{b) at 6q . If we 
denote the difference in the variational ground state energy densities by Aeo = e^^(6|f^) — 
60(60), = fcg + (56 we may write for small g 



Ae,^e^{b^^) + '-^e,{b^^)-eM: 



2 



de^ib) 



db 



m + - 



1 d^e^ib) 



(56)^ + ^^ei(6o^^) (26) 



2 962 
i5br + '-^eM') 



"0 

The first term in the expansion in Sb vanished by definition of the variational minimum. 
Accordingly, the second derivative is positive. The energy ei(6[f^) must be negative to allow 
for (7 > in eq. (24). 



Now we are in the position to decide from the results of the numerical minimization 
procedure whether Acq is bigger or less than zero. We essentially find that g{b^^) = 0{6b), 
typically g = 4- 10-^ ei(6^^) = -O{10-^W), 6b = 5 ■ 10-^ and {d'^e^)/{db^){b^) = 0{W) 
for U/W = 10. Hence, the negative term 5f(6^^)ei(6^^) easily overcomes the positive contri- 
bution (56)^(9^60 )/(962)(6o) by several orders of magnitude. For all U/W we have Acq < 0, 
i.e., the Gutzwiller-Baeriswyl wave function for finite g always wins over the Baeriswyl wave 
function {g = 0). The metal-to-insulator transition again happens at f/?^ = 00. We had to 
do this detailed analysis because the energies of the Baeriswyl and the Gutzwiller-Baeriswyl 
energy become numerically indistinguishable for values of U/W > 5. 

In general, the variational ground state energy density is obtained from a numerical 
minimization of e^^(6, (7) = {T{b, g)) / L + Ud{b,g) with respect to 6 and g. Given the fact 
that the Gutzwiller-Baeriswyl wave function becomes exact both for small and large U/W 
we obtain an excellent estimate for the ground state energy for all U /W. We thus plot 
the relative deviation |(e^^ — eo)/eo| from the exact ground state energy in figure |^. The 
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maximum error is very small, about 1.6% around U/W — 1, and strongly decreases for 
larger and smaller values of the interaction. 

The Gutzwiller-Baeriswyl wave function becomes exact both for small and for large U /W, 
and it gives an excellent estimate for the exact ground state energy density. However, it still 
misses the metal-to- insulator transition in the 1/r- Hubbard model at Uc — W. 

IV. CONCLUSIONS 

In this work we investigated variational approaches to the exactly solvable 1/r-Hubbard 
Hamiltonian. We found a two-parameter trial state, the Gutzwiller-Baeriswyl wave function, 
which becomes exact in both the weak and strong-coupling limit, and gives an excellent upper 
bound for the exact ground state energy. Yet it fails to reproduces the metal-to-insulator 
transition in the half-filled 1/r-Hubbard model at Uc — W. Instead it predicts U^^ — oo, 
just like the Gutzwillcr wave function, = oo. On the other hand, the Baeriswyl wave 
function predicts = 0"*", also in contradiction to the exact result. 

We thus see that even very elaborate correlated wave functions are unable to reproduce 
the correct ground state physics of a one-dimensional model. However, the Gutzwiller- 
Baeriswyl wave function has to be considered an "excellent" wave function according to the 
"energy criterion" . This implies that the "energy criterion" is not very valuable, if one wants 
to draw conclusions about the ground state physics of a given Hamiltonian. This finding is 
supported by the fact that the results of the Gutzwiller Approximation to the Gutzwiller 
wave function are in qualitative agreement with the ground state physics of the 1 / r-Hubbard 
model even though there are considerable quantitative differences. The energy estimate is 
rather poor, and not even variationally controlled, but the basic physical concept that charge- 
correlations drive the transition at a finite interaction strength is properly incorporated in 
the approximation. 

Antiferromagnetic Hartree-Fock theory always overestimates the role of spin-correlations 
in the metal-to-insulator transition while charge-correlations are neglected. Thus, there is 
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no correlation energy for small interactions, the transition takes place for all particle fillings, 
and a finite sublattice magnetization occurs. All three features are not present in the exact 
solution, in part due to the low dimensionality of our system. The physics of the metal- 
to- insulator transition which is predicted to happen at a finite U^^ = W/2 is solely a 
consequence of spin-correlations, contrary to the physics of the transition in the model. 
We know that charge-correlations drive the Mott-Hubbard metal-to-insulator transition, 
and spin-correlations are only residual effects. Thus, the Hartree-Fock picture is always a 
qualitatively wrong description of correlated Fermi systems. 

The situation in higher dimensions is more promising for variational approaches. A 
successful theory of normal fluid ^He is based on the Gutzwiller wave function [^, evalu- 
ated in the limit of high dimensions where the Gutzwiller Approximation to the Gutzwiller 
wave function becomes exact [jl0|,|l3|. In low dimensions, however, systems are Luttinger 
Liquids rather than Fermi liquids. Charge-spin separation in such systems makes it 



difficult to write down qualitatively correct wave functions in terms of the original Fermion 
operators. Since charge-spin separation is common to many one-dimensional systems, it is 
doubtful whether one can write down successful Jastrow-Feenberg wave functions for corre- 
lated itinerant Fermion systems in low dimensions. 

APPENDIX A: CALCULATION OF THE GROUND STATE ENERGY FOR 

CORRELATED WAVE FUNCTIONS 

In ref. we introduced an effective Hamiltonian description of the energy spectrum 
of the 1/r-Hubbard Hamiltonian. We choose a basis in which T is diagonal, and represent 
states by putting hard-core bosons for spin (sjca) and charge {die, e£) degrees of freedom onto 
each site /C, /C = A(m + 1/2) (m = —L/2, . . . , L/2 — 1). They obey a hard core constraint, 
So- ^a: o- + ^k: + ^a: = 1 ^^r each K. Each site is then occupied by one and only one of these 
four objects, t, |, •, o [p^ . 



In this representation the Fermi-sea reads 
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I Fermi-sea) = | [ti] . . . [Ti] 



/c=o 



...o) 



where we assumed an even particle number, and JCp = 7r{2n — 1). To calculate the energy 
expectation value of the Gutzwiller-Baeriswyl wave function we only need to know how T 
and D act on the Fermi sea. From Ref. [^] we know that the kinetic energy operator can be 
split into two parts, 



cS 



-n<K.<TT ^ ^ 



cff 



-7r</C<7r-A 



(Ala) 
(Alb) 



where X/c = sgn(/C), and the prime on the sum in Tq^ indicates that it must not be applied 
to the states represented by [ti] and [•©] which are treated separately by T^^. The total 
effective representation of the kinetic energy operator is T"^^ = Tq^ + T^^. 
The double occupancy operator is expressed by the two terms 



cff 



— 7r</C<7r 
bf= E 



(A2a) 



-7r<K:<7r-A 



i-A^— (2/C + A) 



(A2b) 



— V(27r)2 - (2/C + A)2 4e++^5^,T5^+A,i + H.c. 



The total effective double occupancy operator is D'^^ = t)^ + Df^. 

Fortunately, we only need to apply T'^^ and to the Fermi-sea when we calculate the 
ground state energy of the Gutzwiller-Baeriswyl wave function. The states for /C > ICp are 
already diagonal in T and i), and just give an additive constant. We are left with the action 
of Tf^ and D'^ which only act on neighboring pairs [][]■, and [•o\ on /C, /C + A. We indicate 
this by a prime on the /C-sums. For our purposes we may thus use the following matrix 
representation for 



T 



eff 



T 



K,K+A 



+ E - 



1Cf<IC<tt 



(A3a) 
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fx. ^ 

-Ar 



(A3b) 



and for D'^^ 



D 



cS 



D 



51 Djcjc+A 

-tt<1C<1Cf-A 



\ 



I- at 1 + Xjca/c 



(A4a) 



(A4b) 



where a/c = (2/C + A)/ {2n). The matrix D/c^+a has the eigenvalues zero and one. 

We see that the problem of calculating expectation values (T) and (D) factorizes into 
independent 2x2 matrix problems. We obtain 



ICp<K<iT 



2 
A 



A, 



Ak 



TT<K<K.p-A A/c,/c+aA/c,?C+A 

k:,/c+a-D/c,k:+aA/c,/c+a 



-tt<k<jCf-a A/c,k:+aA/c,/c+a 
Ax;,^+A is the transpose of the vector 



(A5a) 
(A5b) 



A/c,/C+A — Mfcic+A ' 





;:) 








;:) 







for -TT < /C < 



(A6) 



for < /C < /Cf 
with the "correlation matrix" ( 1 is the 2x2 unit matrix) 



V 



(T + ((7-i)z3c,c+A 



b-^>^/^ 

We thus get A/c,/c+a = {a{afc;b,g) (3{aK:]b,g) for -vr < /C < 0, and A^,x;+a 



(AT) 



{(3{aic; b, g) a(a^; b, g) for < IC < ICp where a{x; b, g) and 6, (?) were defined in eqs. (p^G]). 
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Going from /C < to /C > we see that the two entries in A/c,k:+a are exchanged. The 
same happens to row and columns in the Hamilton matrix Hjc^jc+a = T + UD such that 
we may drop the distinction between the two cases. It now is a straightforward task to 
evaluate eqs. ( |A5|) in terms of a{x;b,g) and (3{x;b,g). In the thermodynamical limit one 
easily arrives at eq. (pISD, if one takes into account that the primed sums give an extra factor 
of a half in front of the integral because only every second /C-value has to be included. 

APPENDIX B: ALTERNATIVE CALCULATION OF THE GROUND STATE 
ENERGY FOR THE GUTZWILLER WAVE FUNCTION 

We want to directly calculate the ground state energy of the Gutzwiller wave function by 
using the methods of Vollhardt et al. [|l^,|Tl|]. To this end we first shift the Fermi-sea of filled 
states between k = — tt and kp = Tt{n — 1) by (tt — kp)/2 to generate a symmetric Fermi 
body (rik^cr)- This transformation does not change the double occupancy operator. We may 



then directly apply all formulae of Refs. ||TO||TT|. For example, the formula for the average 



double occupancy d{g) agrees with eq. (|18a|) . To obtain the average kinetic energy we have 



to calculate J2k,a^ [k — {n — kp)/2] (rik^a) where (n^.o-) can recursively be obtained ||T0[- In 
our case we do not need such detailed information to calculate the average kinetic energy. 
We take the lattice periodicity of our dispersion relation into account, and write 

-7r<fc<7r,(T _7r<fc<-7r+(7r-fe^)/2 

where we dropped the term Z)— 7r<fc<7r,cr (uk^a) = which vanishes because the dispersion 
is antisymmetric but the shifted momentum distribution (n^ g.) is symmetric around k = 0. 
The first k-sum in gives the total number of electrons while the second gives half the 
number of electrons outside the Fermi body, nout/2, for which a closed expression can be 
given IjTOl 

1 - c 



'^out — 2- 

^ + 9 



dig) 



(B2) 
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The average kinetic energy then indeed reduces to eq. ( |18b|) . 

We will now show that the Gutzwiller wave function gives the exact ground state energy 
to order {U /W)"^. To this end we recall that the energy spectrum of the kinetic energy 
operator for fixed total momentum consists of equidistant levels, separated by W. The 
potential energy operator can be split into three parts: Dq which commutes with T, and D± 
which scatter eigenstates of T such that the final state has a hj W higher (lower) kinetic 
energy. To first order in U /W we can write the (non yet normalized) exact ground state as 
\iPq) = 1 + at) |Fermi-sea) with a = 0{U /W) because the Fermi-sea is the unique ground 
state (the application of Dq only gives the Fermi-sea again). This expression is nothing but 
the Gutzwiller wave function to order [g — 1) such that the two wave functions agree to first 
order in perturbation theory (the minimization procedure gives 1—g = U /W). Consequently, 
the Gutzwiller wave function gives the exact ground state energy to order {U/WY for the 
1/r-Hubbard model. 

A similar line of arguments can be used to show that the Gutzwiller wave function also 
gives the exact ground state energy to order for all U/W. To this end, we recall that 
the exact ground state wave function for small densities is given by the Fermi sea and 
one particle-hole pair which is created by the action of -D+. For low densities only ladder- 
diagrams contribute which means that no further particle-hole pairs are created but only the 
existing one is scattered. This corresponds to the repeated action of Dq in the subspace of 
the same kinetic energy. This means that all eigenstates of the kinetic energy operator which 
correspond to one particle-hole pair retain the same phase factor. To first order in n we 
may replace all operators Dq, D± by Z). It is thus seen that the exact wave function to first 
order in n has the same states with one particle-hole pair as the Gutzwiller wave function. 
In particular, these states have the same relative phases. Consequently, the Gutzwiller wave 
function is exact to the first non-trivial order in the density, and thus gives the exact ground 
state energy to order for all interaction strength U /W. The corresponding optimal value 
of c/ is g = 1/(1 + U/W). 

We see that lowest-order perturbation theory for the 1/r-Hubbard model is easily done 
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because the Gutzwiller wave function becomes exact in this hmit. 
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FIGURES 

FIG. 1. Hartree-Fock (e^^) and exact (cq) ground state energy density at half-filling as func- 
tion of U/W. For comparison we also include the result of the Gutzwiller Approximation to the 
Gutzwiller wave function, e^"^. 

FIG. 2. Exact (eo) and Gutzwiller /Baeriswyl (cq, Cq ) variational ground state energies as func- 
tion of U/W for the half-filled case; on this scale there is no distinction possible between the exact 
ground state energy density and the variational result from the two-parameter Gutzwiller-Baeriswyl 
wave function, e^^. 

FIG. 3. Relative error Kcg"^ — 6o)/fio| of the Gutzwiller-Baeriswyl variational energy, 6q'^, 
compared to the exact ground state energy density, cq, at half-filling. 
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